A Cellular Automaton Model for the Traffic Flow in Bogota

In this work we propose a car cellular automaton model that reproduces the experimental behavior of traffic flows in Bogot\'a. Our model includes three elements: hysteresis between the acceleration and brake gaps, a delay time in the acceleration, and an instantaneous brake. The parameters of our model were obtained from direct measurements inside a car on motorways in Bogot\'a. Next, we simulated with this model the flux-density fundamental diagram for a single-lane traffic road and compared it with experimental data. Our simulations are in very good agreement with the experimental measurements, not just in the shape of the fundamental diagram, but also in the numerical values for both the road capacity and the density of maximal flux. Our model reproduces, too, the qualitative behavior of shock waves. In addition, our work identifies the periodic boundary conditions as the source of false peaks in the fundamental diagram, when short roads are simulated, that have been also found in previous works. The phase transition between free and congested traffic is also investigated by computing both the relaxation time and the order parameter. Our work shows how different the traffic behavior from one city to another can be, and how important is to determine the model parameters for each city.Comment: 14 pages and 13 figures (gzipped tar file). Submitted to Int.J.Mod.Phys.C. Minor changes, specially at references and typoes, plus a clearer summary of the CA rule

Hodge polynomials of the moduli spaces of pairs

Let $X$ be a smooth projective curve of genus $g\geq 2$ over the complex numbers. A holomorphic pair on $X$ is a couple $(E,\phi)$, where $E$ is a holomorphic bundle over $X$ of rank $n$ and degree $d$, and $\phi\in H^0(E)$ is a holomorphic section. In this paper, we determine the Hodge polynomials of the moduli spaces of rank 2 pairs, using the theory of mixed Hodge structures. We also deal with the case in which $E$ has fixed determinant.Comment: 23 pages, typos added, minor change

Hydrogenated Graphene Nanoribbons for Spintronics

We show how hydrogenation of graphene nanoribbons at small concentrations can open new venues towards carbon-based spintronics applications regardless of any especific edge termination or passivation of the nanoribbons. Density functional theory calculations show that an adsorbed H atom induces a spin density on the surrounding $\pi$ orbitals whose symmetry and degree of localization depends on the distance to the edges of the nanoribbon. As expected for graphene-based systems, these induced magnetic moments interact ferromagnetically or antiferromagnetically depending on the relative adsorption graphene sublattice, but the magnitude of the interactions are found to strongly vary with the position of the H atoms relative to the edges. We also calculate, with the help of the Hubbard model, the transport properties of hydrogenated armchair semiconducting graphene nanoribbons in the diluted regime and show how the exchange coupling between H atoms can be exploited in the design of novel magnetoresistive devices

Molecular dynamics simulations of complex shaped particles using Minkowski operators

The Minkowski operators (addition and substraction of sets in vectorial spaces) has been extensively used for Computer Graphics and Image Processing to represent complex shapes. Here we propose to apply those mathematical concepts to extend the Molecular Dynamics (MD) Methods for simulations with complex-shaped particles. A new concept of Voronoi-Minkowski diagrams is introduced to generate random packings of complex-shaped particles with tunable particle roundness. By extending the classical concept of Verlet list we achieve numerical efficiencies that do not grow quadratically with the body number of sides. Simulations of dissipative granular materials under shear demonstrate that the method complies with the first law of thermodynamics for energy balance.Comment: Submitted to Phys. Rev.

Gyrokinetic and kinetic particle-in-cell simulations of guide-field reconnection. I: Macroscopic effects of the electron flows

In this work, we compare gyrokinetic (GK) and fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields ($b_g$). Here we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma ($\beta_i=0.01$), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field ($b_g\gtrsim 30$) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive ($b_g\gtrsim 5$). Kinetic PIC simulations using guide fields $b_g \lesssim 30$ reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma ($\beta_i=1.0$) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields ($b_g\lesssim 3$).Comment: 18 pages, 13 figures. Revised to match with the published version in Physics of Plasma